Dynamics of saturated Bragg diffraction in a stored light grating in cold atoms 
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We report on a detailed investigation of the dynamics and the saturation of a light grating 
stored in a sample of cold cesium atoms. We employ Bragg diffraction to retrieve the stored optical 
information impressed into the atomic coherence by the incident light fields. The diffracted efficiency 
is studied as a function of the intensities of both writing and reading laser beams. A theoretical model 
is developed to predict the temporal pulse shape of the retrieved signal and compares reasonably 
well with the observed results. 
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I. INTRODUCTION 

The storage of light information in an atomic ensemble 
is a well understood phenomenon which has a promis- 
ing prospect for application both in classical and quan- 
tum information processing Q. The light storage (LS) 
phenomenon allows us to obtain later information on 
a previously stored light pulse, as well as to manipu- 
late the stored information. As it was first proposed, 
LS in an electromagnetically induced transparency (EIT) 
medium is described in terms of a mixed two com- 
ponent light-matter excitation, called dark state polari- 
ton (DSP), where each component of the excitation can 
be externally controlled [3]. To date, several experimen- 
tal observations of these effects were realized in different 
systems 0, i, 0, i, i] ■ 

Alternatively, the LS process can also be described 
as being due to the creation of a spatially dependent 
ground states coherence that contains respectively the 
information on the amplitude and phase of a light pulse 
and which survives after the switching-off of the inci- 
dent light. Using this simpler picture, we have recently 
demonstrated the storage of a polarization light grating 
into an atomic coherence via a backward four-wave mix- 
ing configuration |10|. Others schemes have also been 
recently employed to store spatial structures (images) in 
atomic vapors [ll|, For instance, a light vortex was 
stored in a hot vapor for hundreds of microseconds 

In this work we present experimental and theoretical 
investigation on the dynamics of light grating stored in 
an EIT medium associated with a degenerate two-level 
system. The dependence of the stored light grating with 
the intensities of the incident writing and reading beams 
is investigated. Bragg diffraction into the stored grat- 
ing is employed to probe its dynamics under different 
experimental conditions. The demonstration of the re- 
versible storage and the manipulation of the spatial light 
phase structure stored into the atomic ensemble, and its 
extension to include beams carrying orbital angular mo- 
mentum, would be of great importance to demonstrate 
the manipulation of quantum information encoded in a 
higher dimensional state space [l^ Moreover, the 



storage of this light grating opens up the possibility to 
investigate the generation of correlated photons pairs in 
a previously coherently prepared atomic ensemble jlQ] . 

II. THEORETICAL MODEL 

We consider an ensemble of cold atoms excited by three 
different fields: two writing {W and W') and one reading 
(i?) laser pulses. The atomic ensemble can be well ap- 
proximated by a set of degenerate two-level atoms, with a 
ground-state manifold composed of two degenerate states 
(|la) and |16)) and the excited-state manifold having a 
single state (|2)). As illustrated in Fig. (TJ the ground- 
state degeneracy corresponds to the Zeemam degeneracy 
of atomic cesium in the experiment. In this way, the 
different levels are connected by fields of different polar- 
izations with respect to the atom. We consider fields W 
and R having &^ polarization and field W having ct"^ po- 
larization. W and R excite then the transition lb ^ 2, 
and W the transition la 2. 

The fields W and W propagate in different directions, 
corresponding to a small angle 9 between them. The 
R field is counter-propagating with respect to W. The 
signal we want to model corresponds to the diffraction of 
the R field in the spatial grating formed by fields W and 
W'. In the case of cw excitation of the ensemble, this 
signal corresponds to the well-know conjugated signal in 
four- wave mixing (FWM) processes [lH| . Here we call it 
D field (see Fig.IB)). 

We use this FWM configuration to store and later re- 
trieve a coherence grating written in the atomic ensem- 
ble. In order to address this coherence storage process, 
we use a specific time sequence for the pulsed excitation 
of the ensemble, see Fig.[TJ;. First we prepare the sample 
by exciting it with the two, long writing pulses. In this 
writing process, the goal is to leave the system in its sta- 
tionary state. Then we turn off the writing beams, and 
wait a certain amount of time, the storage time, before 
turning the reading pulse on. This reading pulse stays on 
also for a long time, enough to extract the whole stored 
grating from the ensemble. A field-D pulse is then gener- 
ated during the read process. In the following theoretical 
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analysis, wc want to model and study this field-D gener- 
ation process in detail, considering the three-level-atom 
approximation discussed above. 
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FIG. 1: (a) Simplified Zeeman level scheme, showing the cou- 
pling and the propagation directions of the grating writing 
beams {W and W') and (b) the coupling and the propaga- 
tion direction of the reading (R) and diffracted (D) beams. 
The beams W and W' make a small angle 9 and are cir- 
cularly polarized with opposite handedness, while the beam 
R is counterpropagating to the beam W and have a circu- 
lar polarization opposite to this beam. The diffracted beam 
is detected in a direction opposite to the beam VK'.(c) The 
switching time sequence for the writing and reading beams. 



A. Grating formation and storage 

Consider an atom excited by two writing beams: one, 
Ew, propagating in the z direction and the other, Ew', 
forming an angle 9 with Ew The fields Ew and Ew' 
have orthogonal circular polarizations and (t~ , re- 
spectively. We consider small enough angles so that we 



can assume, to a good approximation, the polarization 
of Ew' as being a~ on the same state basis in which Ew 
is 0-+ . We can then write 



Ew =£w'{r)e"''^'^"^~'^'^'*^cr- 



(la) 
(lb) 



where Ewi^) and £w' {t) represent the transversal modes 
of each field, and we assumed both of them having con- 
stant intensities. The frequencies of the fields arc low 
and uJw'i and their wavevectors are kwz and kw', re- 
spectively. The energy difference between fundamental 
and excited levels is fiLOe. 

The system hamiltonian can then be written as 



H{t) = Ho + V{t) , 

where 

Ho = nLJe\2){2\ 
is the hamiltonian for the free atom and 



(2) 
(3) 



V{t)=-d2,ia-Ew{t)\2){la\ 

-d2,ib-Ew'{t)\2){lb\ + h.c. (4) 

is the interaction hamiltonian. Defining the Rabi fre- 
quencies 



Qw'ir) = 



id2,ia£wir)e 

n 



(5a) 
(5b) 



and assuming the resonance condition ujw = ojw' = '^e, 
the whole set of Bloch equations for the system, in the 
rotating-wave approximation, becomes 



dp2 



dt 

dpia,la 

di 
dpib,ib 
dt 

d<J\a,2 

dt 

daib,2 
dt 

dpia,lb 

dt 



[^w'^ia,2 + ^w'0-ib,2 + C.C.] - T22P22 , (6a) 

[-flwCTla,2 + C.C.] + i„p22 , (6b) 
[-0.W'O-lb,2 + C.C.] + T[f, -^f,p22 , (6c) 

— ^wiP22 — Pla,la) + ^W'Pia,lb — ^120'la,2 , 

(6d) 

— ^W'{P22 — Plb,lb) + ^wPU,la — ^12(7lb,2 , 

(6e) 



= —flw<^2,lb — ^W"7la,2 — 7Plo,16 ) 



(6f) 



with (Tia,2 = /5la.2e-^"»'* and (716,2 = Pife,2 e"*'^*"' * . The 

spontaneous relaxation rates are indicated by and 
r22, for the coherence and population decays, respec- 
tively, r'ljj 1^ and r'j^j indicate the rates at which the 
P22 population decays into the populations pia,ia and 
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Pib,ib, respectively. For simplicity, in these equations and 
in the following, we omit the spatial dependence of the 
Rabi frequencies. The ground-state-coherence decay rate 
7 is introduced to take into account, in an effective way, 
the decay induced by residual magnetic fields. Such de- 
cay is usually a result of inhomogeneous broadening in 
the ensemble of atoms, each subject to a slightly differ- 
ent magnetic field For the signal we are treating 
here, however, this simple model considering the same 
decay constant for the whole ensemble is already enough 
to obtain a good comparison with the experimental data 
probing the coherence decay. 

After a sufficiently long time, the system reaches a 
steady situation in which dpu/dt = 0, for all pu density- 
matrix elements. The steady-state coherence p\g^ be- 
tween the two ground state levels is then given by 



Pla.lb 

with 



^'lb.lb\^w\ 



A 



(7) 



(7r 



12 



■67|17 



(8) 



We are particularly interested in the situation where 7 is 
very small when compared to any other frequency in the 
system, since this corresponds to our experimental con- 
dition. In this limit, note then that the above expression 
simplifies to 



Pla.lb 



w 



(9) 



Once the fields Ew and Ew' are turned off, the coher- 
ences in the system evolve according to their respective 
decay times. Since 7 << ri2, after a time ts >> I/ri2 
the stored coherences in the sample can be well approxi- 
mated by 



<,2(i.)=0, 

<b,2{ts)^0. 
Pla,lh(^s) — PlaAb^ 



(lOa) 
(10b) 
(10c) 



B. Reading 

The stored coherence grating can be extracted from 
the sample using a (t~ -polarized third field En counter- 
propagating with respect to E^y : 

^fl-ffl(f)e^(-'=«^-"«*)<7- , (II) 

with kn, and luji representing the transversal mode, 
wavevector, and frequency, respectively, of field Ef;. Af- 
ter similar considerations as for the grating-formation 



process, including the resonance condition LUe = ujr, and 
the analogous definition of a third Rabi frequency 



(12) 



the relevant Bloch equations describing the reading pro- 
cess become 



daia.2 

dt 

dpiaSb 

dt 



= ^*RPla,lb — ri2Cria,2 , (I3a) 
= ~^Raia,2 - IPla.lb , (13b) 



with aia,2 — Pia,2e~*"^*. Note that the equations for 
o'ia,2 and pia,ib are actually de-coupled from the rest of 
the system of Bloch equations. 

We are interested in calculating the field Ed that is 
phase conjugated to Ew'- This field is generated by the 
medium in the transient excitation of the aia,2 coherence, 
corresponding to the extraction of the stored coherence 
grating. Using the stored state as initial conditions, the 
solution of the above equations for <7ia,2{t) is 

n*RPUibits)e-^''senh{j2t) 

CTla,2(0 = , (14) 
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with 



71 



72 



ri2 + 7 



(ri2-7)'-4|r!i^P 



(I5a) 



(15b) 



The single-atom polarization vector P2,ia on the 2 ^ la 
transition is then given by 



P2,lair,t) = d2.1a(^2,lair,t)e' 



C. Signal 



(16) 



The electric field Ed for the D field coming from the 
diffraction of Er on the sample coherence grating (see 
Fig. [it) is a result of the constructive interference of the 
emission of all atoms in the —kw' direction. If we ne- 
glect interaction between atoms and propagation effects 
on the D field, again for simplicity, the value of Ed in 
the k direction can be obtained by the superposition of 
all atomic contributions on that direction: 



ED{k,t) 



I 



47reo(27r)3/2 



v{r)p2,ia{r.t)e-'''-'^d'r, 



where ry(r) represents the atomic density at r, eo is the 
vacuum permittivity, and the integration runs over the 
whole ensemble volume. Approximating the fields W ^ 
W', and R as plane waves, we can neglect the spatial 
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dependence on £w, 
case, we can write 

ri2 



Ew'i and Er, respectively. In this 



Iw 



21,. 



2Isb 



2Isb 



(18a) 
(18b) 
(18c) 



with Iw, Iw, and Ir the intensities of the W, W' , and R 
fields, respectively. Isa and Ist are the saturation inten- 
sities of the la ^ 2 and 16 ^ 2 transitions, respectively, 
defined according to Ref. [Hj. 

Since ka — kw — 0, Eq. ([TT]) can be written as 



ED{k,t) 



with 



»rf2,la|Pia,lfcl/flWg" 

47reo(27r)3/2 



ri{r)e 



ifci 



(19) 
(20) 



representing the modulus of the stored ground-state co- 
herence, and 



Ir e ^1 *senh (72 t) 



2Is, 



72 /ri 



(21) 



a function describing the temporal profile of the D-field 
pulse. Note that //j(t) is a function of the read field 
parameters only. 

If we approximate the distribution of atoms as having 
a gaussian profile with the same rms width L in all three 
direction, we can write 



r]{r) 



N 



(27rL2)3/2 



-f-r/2L^ 



(22) 



where N is the total number of atoms in the cloud. Using 
this expression for ri{r), Eq. (|19p becomes 



^d2.laN\Pla.lb\fRith' 

47reo(27r)3/2 



'\k+ky 



ED{k,t) 

(23) 

which explicitly shows that the emission of the I?-field 
occurs in the —kw direction only with a spread in vector 
space, on each direction, of the order of the inverse of the 
atomic-distribution spatial width, L"^. 

The detection apparatus can be arranged to collect all 
light in the £)-field mode. In this case, and if the detec- 
tion of the field is performed with a fast detector com- 
pared to the time variation of fR{t), the signal S'fast(^) 
is then proportional to the integration of the intensity of 
light in field D over all k: 



(24) 



5fast(i) =^ / \ED{k,t)\'d'k 



where A is a proportionality constant. From Eq. (|23p 
we see that such detected signal is given by 



(25) 



with A' a different proportionality constant. 

Another important quantity that can be directly de- 
rived from 5'tast(i) is the total energy, Ud, extracted in 
mode D. Note that, in light-storage measurements, the 
goal is usually to extract as much information and en- 
ergy as possible from the coherence grating [l^ . From 
the expressions derived above we have then 



Ud = 



Siast{t)dt 



Ir/2Is. 
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(1 + T^)(2fe + T^) 



(26) 



III. 



EXPERIMENTAL RESULTS AND 
DISCUSSIONS 



As indicated in Fig. la the experiment was performed 
using a degenerate two-level system. This system cor- 
responds in the experiment to the cycling transition 
65i/2(i^ = 3) ^ 6P3/2{F' = 2) of the cesium D2 
line. The cesium atoms were previously cooled in a 
MOT operating in the closed transition 651/2 (^^ = 4) ^ 
6^3/2 (-P"' — 5) with a repumping beam resonant with the 
open transition 65*1/2 (F = 3) ^ 6P3/2(-F' = 3). To pre- 
pare the atoms in the state 65i/2(^ = 3), we switch off 
the repumping beam for a period of about 1 ms to allow 
optical pumping by the trapping beams via non resonant 
excitation to the excited state F' = 4. After optical 
pumping, the optical density of the sample of cold atoms 
in the F — 3 ground state is approximately equal to 3 
for appropriate MOT parameters. 

All the incident laser beams indicated in Fig. 1 are pro- 
vided by an external cavity diode laser which is locked 
to the F = 3 ^ i^' = 2 transition. The grating writing 
beams {W and W) have the same frequency. After pass- 
ing through a pair of acousto-optical modulators (AOM) 
with one of them operating in double passage, they can 
have their frequency scanned around the F — 3 F' = 2 
transition. The two AOM's also allow us to control their 
intensity. These two beams are circularly polarized with 
opposite handedness and are incident in the MOT form- 
ing a small angle 6' w 60 mrad, which leads to a polariza- 
tion grating with a spatial period given by A = 



2sin{e/2) ' 

where A is the light wavelength. The reading beam R 
is circularly polarized opposite to the writing beam W 
and also passes through another pair of AOM's which 
does not change its frequency but allow us to control its 
intensity. 

Employing the time sequence shown in Fig. Ic we have 
investigated the light grating storage dynamics through 
the observation of delayed Bragg diffraction of the read- 
ing beam R in the Zeeman coherence grating induced 
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FIG. 2: Bragg diffraction signal retrieved from the stored 
grating for different storage times 
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FIG. 3: Normalized Bragg diffraction peak signal for different 
storage times. The solid curve corresponds to a fit with an 
exponential function. 



by the writing beams W and W. The writing and 
reading pulses are trigged to the switching off of the 
repumping laser which also triggers the turn off of the 
MOT quadrupole magnetic field. In order to compen- 
sate for spurious magnetic fields, three independent pairs 
of Helmholtz coils with adjustable currents are placed 
around the MOT. 

In Fig. 2 we show the cw-FWM and the Bragg 
diffracted signal which is retrieved from the stored Zee- 
man coherence grating for different storage times. We 
have experimentally verified that the polarization of the 
diffracted beam, both for the steady state cw-FWM sig- 
nal (real time Bragg diffraction) and for the retrieved 
signal, is always opposite to the polarization of the read- 
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FIG. 4: Experimental (right) and calculated (left) retrieved- 
pulsc temporal shape corresponding to different intensities of 
the reading beam, for a fixed storage time. The writing beam 
intensities, W and W' , are 7.0 mW/cm^ and 1.0 mW/cm^ 
respectively. As described in the text, for comparison between 
theory and experiment, all the theoretical reading beam in- 
tensities needs to be re-scaled by a factor a « 0.02. 



ing beam as schematically depicted in Fig. lb. We have 
been able to observe the diffracted signal up to a time 
of 10 /xs. This maximum storage time is very sensitive 
to the compensation of the residual magnetic field. It is 
interesting to note that for short storage times the re- 
trieved signal peak intensity is much larger than the cor- 
responding cw-FWM signal. This effect is related to the 
simultaneous presence of the writing and reading beams 
in the cw regime, where the reading beam contributes to 
decrease the contrast of the coherence grating induced 
by the writing beams. The decay of the peak intensity 
of the diffracted pulse, normalized by its steady state 
value (cw-FWM signal) is presented in Fig. 3. The ex- 
ponential decay behavior is evidenced by the exponential 
fitting (solid curve). For the data presented in Fig. 3, the 
intensities of the writing beams W and W' are approxi- 
mately equal to 5.0 mW/cm^ and 1.5 mW/cm^ respec- 
tively, while the intensity of the reading beam R is about 
8.0 mW/cm^. From the measurement presented in Fig. 
3, we obtain a decay time of the order of 2.9 ^s, which 
corresponds to the Zeeman ground state coherence decay. 
It is worth noticing that we have experimentally verified 
that the measured coherence time does not depend on the 
intensity of either the writing and the reading beams. 

For a fixed storage time of approximately 1/us, we have 
also measured the temporal pulse shape of the retrieved 
signal for different reading beam intensities and the re- 
sults are shown in Fig. 4a-c for three different values 
of the reading beam intensity. We note that the exper- 
imentally retrieved pulse raising time is limited by the 
time constant of the detector (< 0.5 /its). As we have 
discussed previously the coherently prepared atomic sys- 
tem couples to the reading beam to transiently generate 
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the diffracted pulse signal. The temporal width of the 
generated pulse decreases for increasing reading beam in- 
tensity, a direct consequence of the effect of the increased 
dumping of the Zeeman ground state coherence caused by 
spontaneous emission induced by the reading beam in the 
process of mapping the stored Zeeman coherence into the 
optical coherence. In Fig. 4d-f we show the correspond- 
ing retrieved pulse obtained using the previously devel- 
oped theory, assuming Igb the saturation intensity for 
the 6S'i/2(i^ = 3,mF = +3) ^ 6P3/2{P = 2,mF = +2) 
transition. We have used an adjustable parameter of the 
order of a w 0.02 to rc-scalc all the theoretical reading 
beam intensities (i. e., Iji alu), which accounts for 
the uncertainty in the determination of the effective ex- 
perimental value of the Rabi frequency associated with 
the reading beam. 

More systematically, in Fig. 5 we plot the measured 

pulse width (FWHM) for different reading beam intensi- 
ties. In these measurements, for each value of the reading 
beam intensity, we have recorded three curves of the re- 
trieved pulse which allows us to estimate the correspond- 
ing error bars. The solid curve in Fig. 5 corresponds 
to the calculation of the pulse temporal width using the 
signal shape function given by Eq. 25. In this calcu- 
lation we have used 7/ri2 w 0.014 in order to obtain 
the best agreement with the experiment. Note that this 
value is of the same order of the experimentally mea- 
sured decay rate, obtained from the different set of data 
shown in Fig. 3, and estimated as j/ri2 ~ 0.02, with 
2ri2/27r = r22/27r = 5.2 MHz. From the same set of 
data as Fig. 5, we show in Fig. 6 the retrieved pulse en- 
ergy, obtained by time integration of the measured pulse 
intensity. The corresponding solid curve is a theoretical 
fitting obtained using Eq. 26 with the same adjustable 
parameter a. As can be observed the agreement between 
theory and experiment is qualitatively satisfactory own- 
ing to the simplification of the theoretical model, which 
consider a single three-level system and do not accounts 
for the manyfold Zeeman degeneracy. In fact, for the 
intensities of the writing beams used in the experiment, 
ground state coherence involving different pairs of Zee- 
man sub-levels can actually exist. 

We also have measured the variation of the diffracted 
signal as a function of the intensity of one of the grat- 
ing writing beam (i.e., the beam W) and the results for 
the corresponding pulse energy are shown in Fig. 7. For 
these measurements, the intensities of the grating writing 
beam W' and the reading beam were respectively equals 
to 1.0 mW/cm^ and 9.0 mW/cm^. The solid curve in 
Fig. 7 corresponds to a theoretical fitting with the cal- 
culated retrieved pulse energy given by Eq. 26, assum- 
ing that Isa is the saturation intensity of the 6S'i/2(-F' = 
3,mF = +1) — * 6P3/2(-F' = 2,rnF = +2) transition, with 
Isa = 15/s6 according to the ratio between the corre- 
sponding Clebsch-Gordan coefficient. Again, to account 
for the uncertainty in the experimental value of the Rabi 
frequency associated with the writing beams W and W , 
we have used another adjustable parameter, which in 
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FIG. 5: Measurement of the temporal width (Full Width at 
Half Maximum) of the retrieved pulse for different intensities 
of the reading beam, obtained in similar experimental condi- 
tions as in Fig. 4a-c. The solid curve is a theoretical fitting 
using the model described in the text. 
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FIG. 6: Retrieved pulse energy for different intensities of the 
reading beam, obtained on similar experimental conditions 
as in Fig. 4(a-c). Inset: The corresponding variation of the 
peak intensity of the retrieved pulse. The solid curves are 
theoretical fittings employing the model described in the text 
with the same intensity adjust parameter used to fit the pulse 
width in Fig. 5. 

the present case is of the order of a! k, 1.9, to re-scale 
the theoretical intensity ratio between these beams (i.e., 

Iw jlw a' hv / Iw) ■ 

As can be observed from Fig. 6 and Fig. 7 the amount 
of energy that can be retrieved from the medium clearly 
saturates with the writing and reading beam intensities. 
In particular, this shows that for fixed writing beams 
intensities, there is a maximum amount of energy that 
can be retrieved from the stored coherence. However, 
it is worth mentioning the different saturation behavior 
observed for the variation of the peak intensity of the re- 
trieved pulse with the corresponding intensity, as shown 
in the insets of Fig. 6 and Fig. 7. As observed, the maxi- 
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Writing beam intensity, W (mW/cm ) 

FIG. 7: Variation of the maximum peak intensity of the re- 
trieved pulse for different intensities of the grating writing 
beam W. For these data, the corresponding intensities of 
the writing (W) and reading (R) beams were fixed at 1.0 
mW/cm^ and 9.0 mW/cm^, respectively. The solid curves 
are again theoretical fittings using the model described in 
the text. We have used the same intensity adjust parame- 
ter, a' « 1.9, in both curves. 

mum peak of the retrieved pulse saturates more strongly 
with the writing beam intensity as compared with the 
saturation induced by the reading beam. As one should 
expect, the saturation induced by the reading beam is 
related mainly to the total retrieved energy. The pulse 
peak, however, can increase much further with the read- 
ing power, since it is closely related also with the speed 
of the reading process. On the other hand, the increase 
of the writing beam intensity will saturate the Zeeman 
coherence grating, therefore reducing its contrast. This 



effect has a strong influence on the Bragg diffraction ef- 
ficiency, affecting equally the total retrieved energy and 
the pulse peak. 



IV. SUMMARY 

We have investigated, both theoretical and experimen- 
tally, the storage of a spatial light polarization grating 
into the Zeeman ground states coherence of cold cesium 
atoms. Systematic measurements were performed to re- 
veal the saturation behavior of the retrieved signal as 
a function of the intensities of the writing and reading 
beams. The developed simple theoretical model accounts 
reasonably well for the observed results and in partic- 
ular for the measured pulse temporal shape. We con- 
sider our results are of considerable importance for a 
better understanding of the coherent memory for mul- 
tidimensional state spaces. Finally, we would like to 
mention that we also have observed the coherent evo- 
lution of the stored grating in the presence of an applied 
magnetic field, which shows collapses and the revivals 
of the stored coherence grating. This effect is associated 
with the Larnior precession of the induced grating around 
the app lied magnetic field as was reported previously in 
[20I . l2l( and strongly support the possibility of manipu- 
lating more complex spatial information stored into an 
atomic medium. Further investigation on this effect is 
currently under way and will be published elsewhere. 
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